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Wc present a measurement of the mass of the top quark using data corresponding to an integrated 
luminosity of 1.9 fb~^ of pp collisions collected at ^/s = 1.96 TeV with the CDF II detector at 
Fermilab's Tevatron. This is the first measurement of the top quark mass using top-antitop pair 
candidate events in the lepton + jets and dilepton decay channels simultaneously. We reconstruct 
two observables in each channel and use a non-parametric kernel density estimation technique to 
derive two-dimensional probability density functions from simulated signal and background samples. 
The observables are the top quark mass and the invariant mass of two jets from the W decay in 
the lepton + jets channel, and the top quark mass and the scalar sum of transverse energy of the 
event in the dilepton channel. We perform a simultaneous fit for the top quark mass and the jet 
energy scale, which is constrained in situ by the hadronic W boson mass. Using 332 lepton + 
jets candidate events and 144 dilepton candidate events, we measure the top quark mass to be 
Mtop = 171.9 ± 1.7 (stat. + JES) ± 1.1 (other syst.) GeV/c^ = 171.9 ± 2.0 GeV/c^ 
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I. INTRODUCTION 

With a mass of approximately 172 GeV/c^ [T], the top 
quark (t) is by far the most massive fundamental object 
observed to date in nature, some 40 times as massive as 
its isospin partner, the bottom quark (6). This large mass 
leads to an important role for the top quark in theoreti- 
cal predictions from the standard model (SM) of particle 
physics. In particular, electroweak radiative corrections 
to the W boson mass, due to loops containing top quarks, 
play an important role in constraining the mass of the 
Higgs boson, which also contributes to radiative correc- 
tions. If the Higgs boson is discovered, a precise mea- 
surement of the mass of the top quark will help provide 
an important test of the SM, and would confirm that the 
newly observed object is the SM Higgs boson and not 
some other scalar particle or source of new physics. In- 
dependent of the Higgs boson, the large mass of the top 
quark may make precision measurements throughout the 
top quark sector necessary to help disentangle models of 
new physics 0. 

The CDF and D0 collaborations jointly announced 
discovery of the top quark in 1995 [1, but it was 
not until the availability of large datasets from Run 
II at the Tevatron that precision measurements of the 
top quark mass have been possible. The 2007 world 
average of published top quark mass measurements, 
Mtop = 172.5 ± 1.5 (stat.) ± 2.3 (syst.) GeV/c^ [l|, 
compares to the CDF measurement of Mtop — 174 ± 
10 (stat.) t\l (systj GeV/c^ upon finding first evidence 
for the top quark According to the SM, top quarks 
at the Tevatron are produced mainly in tt pairs result- 
ing from qq annihilation (85%) and gluon-gluon fusion 
(15%). A top quark decays more than 99% of the time 
to a, W boson and a b quark. The topology of tt events 
depends on the subsequent decay of each of the two W 
bosons. Each W boson can decay hadronically, to a pair 
of quarks, or leptonically, to a charged lepton and a neu- 
trino. Due to the difficulty of reconstructing r leptons 
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when they decay, for the purposes of this analysis only 
electrons (e) and muons (/x) are considered as charged 
lepton candidates. The semileptonically decaying r lep- 
tons can enter the data set when they are reconstructed 
as electrons or muons. 

Dilepton tt events are those in which both W bosons 
decay leptonically. The dilepton channel has a small 
branching ratio of 5%, and suffers from undercon- 
strained kinematics resulting from the presence of two 
neutrinos in each event that escape undetected. The ad- 
vantages of the dilepton channel are a low background 
rate and simple combinatorics with only two quarks in 
the final state. Lepton + jets events are the roughly 
30% of tt events in which one W boson decays hadroni- 
cally and the other decays leptonically. The lepton -I- jets 
channel has four quarks, one lepton, and one neutrino in 
the final state, and sufficient amount of information mea- 
sured in the detector to constrain the kinematics of the 
tt decay. 

This analysis describes a measurement of the top quark 
mass in both the lepton -I- jets and dilepton decay chan- 
nels using data collected at the Tevatron with the CDF 
II detector in a 1.9 fb^^ integrated luminosity nm. This 
is the first analysis to combine likelihoods from multi- 
ple measurements of the top quark mass in different de- 
cay topologies into a single joint likelihood with a robust 
treatment of the correlations in systematic uncertainties 
between the two channels channels. 

In the analysis described in this article we follow the 
template strategy Q. In the lepton -I- jets channel we 
determine the kinematics of the decay by fitting for a 
reconstructed top quark mass ml'^'^° In the dilep- 
ton channel, due to the two undetected neutrinos, there 
is not enough information to constrain the four-vectors 
of the top quarks. Instead we use the neutrino weight- 
ing algorithm (NWA), in which we scan over top quark 
masses, performing an integration over the polar angles 
of the neutrinos at each mass to calculate a weight based 
on the agreement with the measured momentum imbal- 
ance in the event [1, [§] . We select as the observable in an 
event the top quark mass (mf^-^) that yields the highest 
weight. 

Uncertainties in jet modeling and in the calorimeter 
response result in a systematic uncertainty in the jet en- 
ergy calibration, which in turn induces the largest sys- 
tematic uncertainty in top quark mass measurements. In 
the lepton -I- jets channel, the hadronically decaying W 
boson provides an in situ calibration sample for these ef- 
fects. The invariant mass of the hadronically decaying 
W bosons {iTijj) is used in the likelihood fit to measure 
and constrain the jet energy calibration. This procedure 
reduces the combined statistical and systematic uncer- 
tainty on the top quark mass measurement. The cal- 
ibration obtained from the lepton -I- jets channel is for 
the first time applied to both channels by performing a si- 
multaneous likelihood fit. To improve the precision of the 
dilepton measurement we also use a second observable, 
Ht which is a scalar sum of the ti decay product 
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transverse energies. 

Conventional procedures for combining values of the 
top quark mass measured in different decay channels re- 
quire correlations in the systematic effects between the 
channels as inputs [T]| . In addition, combinations of ex- 
isting top quark mass results must also assume the func- 
tional form of the likelihood in each measurement. Typ- 
ically, the likelihoods are assumed to be Gaussian, and 
likelihood asymmetries or other departures from Gaus- 
sian behavior are not taken into account. The analysis we 
describe here is free of these two assumptions. The sys- 
tematic uncertainties are evaluated using the simultane- 
ous fit, and the measurement uses one likelihood function 
that depends on the data in the two channels. As a cross- 
check, we also measure the top quark mass separately in 
the lepton -I- jets and dilcpton channels, including a full 
evaluation of systematic uncertainties. 



II. EXPERIMENTAL SETUP 
A. CDF II Detector 

The Collider Detector at Fermilab is located at one of 
two collision points along the ring of the Tevatron acceler- 
ator, which collides bunches of protons and anti-protons 
at a center-of-mass energy of 1.96 TeV. The detector has 
an approximate cylindrical geometry around the Teva- 
tron beamline and is described in a Cartesian or in a 
polar coordinate system. In Cartesian coordinates the z 
axis is located along the beam axis with positive z in the 
direction of the proton beam, the x axis pointing outward 
in the plane of the Tevatron ring, and the y coordinate 
pointing up. It is often more convenient to use polar co- 
ordinates: the azimuthal angle is the angle from the 
X axis in the plane transverse to the beamline; the po- 
lar angle Q is the angle from the proton beam direction. 
The pseudorapidity rj = — ln(tan |) is a quantity numer- 
ically close to rapidity for highly relativistic particles; 
differences in pseudorapidity are therefore nearly invari- 
ant with respect to boosts along the z axis. Collisions 
occur along the beamline and are distributed about the 
center of detector with a spread of about 30 cm. We dis- 
tinguish between ij defined with respect to z = (rjdet) 
and rj defined with respect to the event collision point. It 
is common to reference a region in a cone of AR around 
an object. This refers to the nearby region in rj-cj) space: 
AR = \J {Aj])'^ + (A0)2. The transverse momentum p-r 
refers to the momentum in the plane transverse to the 
beamline. The transverse energy of an object is defined 
as Et = EsmQ. 

A detailed description of the CDF II detector is pro- 
vided in Ref. [12]. In this section we briefly introduce 
the detector subsystems relevant to this analysis, start- 
ing with the detectors closest to the interaction region. 

Charged particles are observed in the silicon tracking 
detectors. The innermost silicon detector, layer 00 is 
a single-sided silicon strip detector mounted directly onto 



the beryllium beampipe, providing axial tracking infor- 
mation at a radius of 1.6 cm. The silicon vertex detector 
(SVX II) consists of five double-sided silicon strip detec- 
tors locating at radii from 2.5 cm to 10.6 cm from the 
beamline and 90 cm in length, providing axial and stereo 
information. Tracking of charged particles in the central 
region (iTydetl < 1-0) is provided by a 310 cm long cylin- 
drical open cell drift chamber, the central outer tracker 
(COT), located at radii between 43 and 132 cm. The 
tracking detectors are immersed in a 1.4 T solenoidal 
magnetic field, allowing for charge determination and 
momentum measurements of charged particles p^ . 

The calorimeter system measures the energy and po- 
sition of particles passing through and interacting with 
dense material. CDF uses lead-scintillator and steel- 
scintillator sampling devices for the electromagnetic 
and hadronic calorimetry, respectively. The calorim- 
eter system is comprised of the central electromagnetic 
(CEM) central (CHA) and wall (WHA) hadronic 

calorimeters [l^ covering \T\det\ < 1.0, and the plug elec- 
tromagnetic (PEM) [TtI and hadronic (PHA) calorim- 
eters covering 1.1 < \r\det\ < 3.6. Shower maximum 
detectors are embedded in the central (CES) and plug 
(PES) electromagnetic calorimeters at approximately six 
radiation lengths from the collision point to provide the 
transverse shape of the shower jT^, [ll| . 

Muon detectors are located beyond the calorimeters. 
Directly outside the CHA is the central muon detector 
(CMU), which covers \r\det\ < 0.6. Located behind fur- 
ther 60 cm of steel shielding is the central muon upgrade 
(CMP) detector. The central muon extension (CMX) 
covers the region 0.6 < |?7det| < 1-0. Muons pass through 
the calorimeter and the shielding and leave behind a se- 
ries of hits (stubs) in the muon detectors, which consist 
of four layers of single- wire drift cells. 

CDF employs a three-level trigger system to select 
potentially interesting events, reducing the interaction 
rate from the 1.7 MHz average bunch crossing rate to a 
more manageable 75-150 Hz. The tt candidate events 
used in this analysis are collected by triggers that iden- 
tify at least one high-py lepton candidate. For cen- 
tral electron events (CEM events), the first-level trig- 
ger requires a cluster in the electromagnetic calorimeter 
with Et > 8 GeV, a matching track in the COT with 
Pt > 8 GeV/c, and a ratio of energy deposited in the 
hadronic to electromagnetic calorimeters less than 1:8. 
At the second trigger level, the cluster energy require- 
ment is tightened to Et > 16 GeV, and the third trigger 
level makes basic electron identification cuts and further 
tightens the energy requirement to Et > 18 GeV. For 
muons in the central region, the first-level trigger requires 
stubs in both the CMU and CMP detectors (CMUP 
events) or the CMX detector (CMX events) and a match- 
ing track in the drift chamber with pT > 4 GeV/c for 
CMUP events and with pT > 8 GeV/c for CMX events. 
At the second trigger level, the pt requirement is in- 
creased to 15 GeV/c, and at the third trigger level to 
18 GeV/c. 
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III. EVENT SELECTION AND BACKGROUND 
ESTIMATION 

The energetic, charged leptons and missing transverse 
energy from at least one leptonic W boson decay help 
distinguish lepton + jets and dilepton tt events from the 
QCD multijet background. Further rejection of back- 
ground events containing real leptons is achieved by re- 
quiring high-pT jets, and in some cases identifying one or 
more of those jets as arising from a b quark. We briefly 
describe here the reconstruction of physics objects in the 
detector, as well as event selection, background estima- 
tion, and Monte Carlo (MC) simulation. 



A. Selection of leptons 

Lepton identification is similar in the lepton + jets 
and dilepton channels. Events in both channels require 
a clean lepton in the central region of the detector. The 
second lepton in dilepton events can be from less pure 
categories such as forward or non-isolated leptons. The 
major distinctions between the two channels are noted in 
the following descriptions. 

A small cluster of towers in the CEM containing Et > 
20 GeV with a COT track that extrapolates to the face 
of a tower in the cluster is identified as an electron candi- 
date. Electrons deposit most of their energy in the elec- 
tromagnetic calorimeter; therefore we require the ratio 
of hadronic energy to electromagnetic energy to be less 
than 0.055-|-0.00045i?, where E is the total energy of the 
electron. To reject backgrounds with energetic tt^s plus a 
track, we require the ratio of energy in the cluster to the 
track momentum to be not more than 2.0 for electrons 
with Et < 100 GeV. The lateral shower development 
measured in the calorimeter and the CES is required to 
match the electron shower shape as measured on a test 
beam. We also require that the COT track extrapolated 
to the depth of the shower maximum detector matches a 
CES cluster in the r-z and r-(p planes p^ . 

Full COT tracking information is available only for 
IVdetl < 1-0. To reconstruct forward (PHX) electron can- 
didates with |?7det| > 1.0, we use clusters in the PES, 
the energy measurement in the PEM, and knowledge of 
the pp interaction vertex. Silicon detector hits are then 
added to form a track . We require that forward elec- 
trons have Et > 20 GeV and a ratio of hadronic to elec- 
tromagnetic energy less than 0.05. We also require that 
shower profiles measured in the PEM and PES match the 
electron shower shape as measured on a test beam 19]. 
While PHX electrons are accepted in the dilepton sam- 
ple, they are used in the lepton -I- jet sample only to veto 
dilepton events. 

Muon candidates are required to have a track with 
Pt > 20 GeV/c matched to a calorimeter tower with elec- 
tromagnetic energy less than 2 GeV and hadronic energy 
less than 6 GeV. The energy of the tower is required to 
be low since muons are minimum ionizing in the calo- 



rimeter material. Both requirements are slightly looser 
if the track momentum is greater than 100 GeV/c. If 
the track extrapolates to a fiducial region of the muon 
chambers, we require that it is matched to a stub inside 
these detectors. We categorize muons based on the de- 
tector chamber that the muon traverses. Muons fiducial 
to both the CMU and CMP detectors are called CMUP 
muons. If a track does not extrapolate to any muon 
chambers but all other quality criteria are satisfied, we 
accept this object as a "stubless" muon (CMIO). Only 
CMUP and CMX muons enter the lepton -I- jets dataset. 
All muon categories are allowed into the dilepton dataset 
for one of the two leptons to increase the statistics of the 
sample. 

Electrons and muons produced in W boson decays will 
in general be well separated from other objects in the 
event. For electrons we define isolation as the ratio of 
energy deposited within AR < 0.4 around and excluding 
the electron cluster to the electron cluster energy itself. 
For muons, isolation is defined as the ratio of the trans- 
verse energy within AR < 0.4 from the tower crossed by 
the muon to the muon pT- The tower traversed by the 
muon is excluded from transverse energy sum. We cate- 
gorize leptons as isolated if their isolation variable is less 
than 0.1. The lepton -I- jets channel uses only isolated 
leptons. In the dilepton channel, to increase the sample 
statistics we allow one non-isolated lepton of any type 
except for PHX electrons and CMIO muons. 



B. Jet corrections and systematic uncertainties 

The property of quark confinement [2l|, [2^ ensures 
that bare quarks are not directly observable after QCD 
processes take effect. Quarks and gluons ("partons") in- 
stead manifest themselves in the detector as jets of parti- 
cles flowing in the direction of the original gluon or quark 
(parton). The sum of energies of particles within a cone 
around the direction of the fragmenting parton is strongly 
correlated with its energy. Jets for this analysis are 
reconstructed with the cone-based clustering algorithm 
JETCLU [2^ , using a cone in the azimuth-pseudorapidity 
space AR = y'rf + cj)^ = 0.4. The four-vector of a jet 
is constructed based on energies and locations of calo- 
rimeter towers belonging to it. A detailed explanation 
of CDF calibration of jet energy, the corrections applied 
to it and the associated systematic uncertainties can be 
found in Ref. We briefly summarize these correc- 

tions below. 

From the raw jet energy several stages of corrections 
are applied that attempt to remove effects masking the 
initial parton energy. In the first stage, the relative cor- 
rections normalize the detector response as a function of 
i]det so that jets at all r]det have the same Et response as 
jets in the well understood central region of the detector, 
0.2 < |77det| < 0.6. The response varies across ry^et due 
to uninstrumented regions of the detector and different 
amounts of material in front of the calorimeters. The 
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FIG. 1: Fractional uncertainty on corrected jet pt as a func- 
tion of corrected jet pr- 



energy (^t) is defined with respect to the event's inter- 
action vertex: 



-^7 



(1) 



where the sum on i runs over all calorimeter towers, 
is a unit vector in the transverse plane pointing from 
the beamline to the face of the ith calorimeter tower, 
and is the Et in the ith tower, calculated using the 
reconstructed event vertex. In events with muons, the 
measured pt of the muon track is used in the calcu- 
lation, and not the energy deposited by the muon in the 
calorimeters. Since jets are measured with better energy 
resolution than the raw calorimeter towers, calorimeter 
tower energies are replaced with jet energies for the tow- 
ers clustered into jets. 



next correction accounts for the average energy in a jet 
cone due to additional pp inelastic interactions occurring 
during the same bunch crossing. After accounting for 
energy due to extra pp interactions, corrections are ap- 
plied for calorimeter non-linear response to hadron jets 
and energy loss in uninstrumented regions. At this stage, 
jet energies should be independent of the CDF detector, 
and correspond to the energies of "particle jets" , which 
are defined as all long lived particles from the primary 
pp collision within the jet cone. Energy from spectator 
partons in the hard collision process that breaks up the 
proton and anti-proton to form the tt system is accounted 
for by the underlying energy correction. The final cor- 
rection accounts for out-of-cone effects, in which some of 
the original parton energy lies outside the jet cone. 

Modeling of each of the effects described above is a 
potential source of uncertainty on the measurement of 
jet energies. The combined fractional uncertainty on the 
jet energy calibration (ctc) is shown in Fig.[TJ The overall 
jet energy calibration is referred to as the jet energy scale 
( JES) . We measure the difference A jes between the JES 
effects in simulation and data in units of ctc. The a priori 
AjEs estimate at CDF is thus by definition 0± 1 CTc. For 
the lepton -|- jets channel, jets with \ridet \ < 2.0 and Et > 
20 GeV after applying the relative, multiple interactions, 
and hadron jet response linearity corrections are referred 
to as "tight jets." Jets not passing the tight cuts but 
having Et > 12 GeV and |?ydet| < 2.4 are referred to as 
"loose jets." The dilepton channel uses jets with \r]det \ < 
2.5 and Et > 15 GeV after corrections to the particle jet 
level. 



C. Neutrinos 

Neutrinos from leptonic W decays escape the detec- 
tor undetected, leading to an imbalance of energy in the 
transverse plane of the detector. This missing transverse 



D. Identification of b jets 

The B hadrons produced in the fragmentation of b 
quarks have an average lifetime on the order of 1.5 ps. 
Given typical boosts of b jets in tt events at the Teva- 
tron, this implies that B hadrons on average travel sev- 
eral millimeters in the detector before decaying, leading 
to the identifiable signature of a displaced vertex. This 
analysis uses the SECVTX algorithm [2^ to find jets con- 
taining secondary vertices consistent with the decay of a 
B hadron; when such a displaced vertex is found, the jet 
is identified ("tagged") as a 6-jet. In the lepton -I- jets 
(dilepton) channel, only the four (two) jets with largest 
Et are checked for b tags. Only tight jets are allowed 
to have 6 tags in the lepton -|- jets channel. The tagging 
efficiency of b quark jet depends on the Et and rjdet of 
the jet, but is roughly 40% within the fiducial acceptance 
of the silicon detector, with a roughly 1% per-jet prob- 
ability to incorrectly tag a light flavor jet. Jets arising 
from charm quarks have finite lifetimes as well, and have 
a per-jet tagging efficiency of roughly 8%. 

E. Lepton + jets selection and background 
estimation 

The lepton -|- jets decay channel has four quarks in the 
final state, so we require at least four jets in the detector. 
Only events with at least one b tag are considered; events 
with no b tags have a substantial background from pro- 
duction of W bosons with four or more jets, and are also 
particularly sensitive to systematic uncertainties. We di- 
vide lepton -|- jets events into two exclusive subsamples 
based on the number of b tags. For events with a single 
b tag (1-tag), exactly four tight jets are required. Events 
with additional tight jets and events with only three tight 
jets, but one or more loose jets were predicted by MC sim- 
ulation to contribute very little in terms of measurement 
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TABLE I: Lepton + jets event selection summary. Jets are 
corrected to the particle jet level. 



1-tag 2-tag 

6-tags (Leading 4 jets) =1 > 1 

Lepton pT (GeV/c), Et (GeV) > 20 

(GeV) > 20 

Leading 3 jets Et (GeV) > 20 

4tli jet St (GeV) >20 >12 

Extra jets Et (GeV) <20 Any 



precision. Events with two or more tags (2-tag) have ex- 
tremely small background contamination, and also have 
fewer ways to assign jets to the quarks from the tt de- 
cay, so we allow events with three tight jets and one or 
more loose jets, as well as events with four or more tight 
jets. For each lepton -I- jets subsample, we also require 
a single, high-pr, isolated CEM, CMUP or CMX lepton 
and higher than 20 GeV, indicating the presence of 
an escaping neutrino. Events with multiple leptons are 
vetoed. The lepton -|- jets event selection is summarized 
in Table m 

Background estimates for the lepton -I- jets events are 
derived from a hybrid of data- and MC-based measure- 
ments, similarly to previous dedicated analysis [2^. Data 
are used whenever possible, and MC information is used 
to fill in any remaining gaps in knowledge. In particular, 
the overall rate of events with real W bosons and addi- 
tional jets {W -I- jets), which dominate the background 
sample, is determined using the data. The fractions of 
these events with one or two charm quarks ( Wc and 
Wcc events) and two bottom quarks {Whb events) are 
determined from MC samples. Overall normalizations 
of -I- jets events come from the data after subtract- 
ing off an estimate for the fraction of W events coming 
from QCD processes (non-VF events), and separating out 
a MC based estimate for other processes with real W^'s 
{WW/WZ/ZZ, tt and single-top production), where the 
lepton trigger efficiencies and MC-data differences have 
been taken into account. The number of background 
events also depends on the rate to mistag light-quark 
jets; this fake rate is determined using data samples trig- 
gered by presence of jets. Fake tags come mostly from 
mismeasured tracks, through interactions with material 
in the detector and real decays of long-lived light-flavor 
particles such as Kg and A also contribute. 

Energetic charged leptons from W boson decay can 
be faked by QCD events via conversions (electrons) or 
misidentified pions and kaons (muons), as well as from 
semi-leptonic heavy-flavor decays. In such events, the 
-^T requirement can also be passed when jets are mismea- 
sured or fall into uninstrumented regions of the detector. 
The events which do not include a W boson, typically do 
fail the requirement, however, so the entire dis- 
tribution, including the Iow-^t region, is used to fit for 
the number of events with fake W bosons and high ^t- 



TABLE II: Expected event yield for the lepton -I- jets selection 
after all cuts. Uncertainties quoted are due to the uncertainty 
on the integrated luminosity, statistics of MC samples and 
uncertainty in QCD modeling. 





1-tag 


2-tag 


Wbb 


9.1±3.7 


2.1±0.9 


Wcc, Wc 


8.3±3.4 


0.5±0.3 


W (mistags) 


10.4±2.3 


0.2±0.1 


Single top 


2.0±0.1 


0.7±0.1 


Diboson 


2.4±0.2 


0.21±0.02 


QCD 


10.4±8.7 


0.3±1.6 


Total Background 


42.7 ± 12.5 


4.2 ± 1.9 


tt (6.7 pb) 


156.7±21.1 


76.6±12.0 


Observed 


233 


99 



The QCD background is modeled in the data by events 
passing all cuts where instead of a lepton an electron-like 
object is required. These pass all the kinematic cuts that 
are imposed on the electrons, but fail shower develop- 
ment or track quality cuts. In an alternative, high statis- 
tics model for the QCD background data events are used 
where the isolation cut on the charged lepton is removed 
and an anti-isolation cut is imposed, requiring the isola- 
tion variable defined earlier to be greater than 0.2. The 
non-isolated model is not used to obtain the QCD back- 
ground normalization, however it is used at later stages 
in the analysis. 

Table [iTl shows the expected number of background 
events in the 1-tag and 2-tag samples after all cuts, the 
expected number of signal events based on the theoret- 
ical cross-section (6.7 pb) ^ at Mtop = 175 GeV/c^, 
and the observed number of events. Uncertainties on the 
event yield expectations are due to integrated luminos- 
ity, statistics of the MC samples as well as uncertain 
spectrum in the QCD model. 

F. Dilepton selection and background estimation 

Signal events in the dilepton channel contain two b 
quarks, two neutrinos and two oppositely charged lep- 
tons. We require that at least one of the leptons is iso- 
lated and categorized as a CEM, CMUP or CMX lepton. 
The second lepton can belong to any category. If the 
second lepton is a PHX electron or a CMIO muon, it 
must be isolated. All other types of second leptons can 
be non-isolated to increase statistics of the sample. We 
require a minimum of two jets in each dilepton event. 

To account for the two neutrinos, we require > 
25 GeV. We increase the requirement to > 50 GeV 
if the angle in the r-(j> plane between the vector and 
any jet or lepton is less than 20° in order to reduce back- 
grounds such as Drell-Yan production of r pairs and QCD 
events where jets fall into uninstrumented regions of the 
calorimeter. 

We require that the Ht - the scalar sum of transverse 
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energies of jets, leptons and be greater than 200 GeV. 
This has a small effect on signal acceptance, as the sum 
of energies of the tt decay products must be equal to or 
greater than twice the top quark mass. 

In events where the charged leptons are the same fla- 
vor, backgrounds with Z bosons are removed by requiring 
the invariant mass of the dilepton pair to be smaller than 
76 GeV/c^ or larger than 106 GeV/c^. We impose this 
requirement only if the significance (s^^) is smaller 
than 4.0GeV^/^: 



sum, 
T 



< 



4.0 GeV^/2 



(2) 



where is defined as the scalar sum of the transverse 

energies of all calorimeter towers, with muons and jets 
corrected as in the calculation. 

Backgrounds for the dilepton channel include DrcU- 
Yan processes, diboson production and QCD multijet 
production. We estimate the contribution from Drell- 
Yan production of r pairs and diboson production using 
MC samples normalized to the theoretical cross sections. 
We apply trigger efficiencies as well as corrections ac- 
counting for differences between the data and MC simu- 
lations in lepton identification efficiencies and jet multi- 
plicity distributions. 

We employ a combined data-MC sample technique to 
estimate the contamination from Drell-Yan ee and fifi. 
This background contains two components in the sig- 
nal region: events outside the Z boson window (76 — 
106 GeV/c^), and events inside the window and pass- 
ing the significance cut. We first count the number 
of events in data with a dilepton invariant mass within 
the Z window. After subtracting expected contributions 
from other sources, we multiply this estimate by the ratio 
of the number of events outside the window to the num- 
ber of events inside the window, as measured in the MC 
samples. This gives the number of Drell-Yan events out- 
side the Z window. We estimate the contribution from 
events inside the Z window with high by multiply- 
ing the number of events in data inside the window by 
the ratio of events passing and failing the M^t significance 
cut, again obtained in MC samples. 

The data are used to estimate the contribution of 
events where a real lepton is produced in association 
with multiple jets and one of the jets is misidentified as 
a second lepton. Data samples triggered on presence of 
jets are used to obtain the probability for a jet to fake 
a charged lepton. These probabilities depend on lepton 
category and the jet Et- We apply these probabilities to 
the single-lepton data to obtain an estimate for the fake 
background contribution. 

We divide the dilepton sample into non-tagged and 
tagged subsamples, which have very different purity. 
Since the fake background is modeled directly from data, 
the probabilities for fake leptons to be reconstructed are 
summed separately in events with and without a b tag. 
All other backgrounds are modeled using MC samples. 
For all MC events, we calculate the probability for each 



TABLE III: Expected event yield for the dilepton selection 
after all cuts. Uncertainties quoted capture the uncertainty 
on integrated luminosity, statistics of the MC samples and 
uncertainties on the fake rates. 





non-tagged 


tagged 


Diboson 


9.1 ±2.2 


0.3 ±0.1 


Drell-Yan 


16.0 ±2.5 


0.9 ±0.1 


Fakes 


19.3 ±5.6 


2.7 ±1.0 


Total Background 


44.3 ±7.0 


3.9 ±1.0 


tt (6.7 pb) 


40.1 ±3.1 


55.8 ±4.2 


Observed 


83 


61 



jet to be tagged, accounting for the probabilities for light 
flavor jets to be mistagged. Given the tag probabilities 
for the two leading jets, we calculate the probability for 
each event to enter the non-tagged and tagged subsam- 
ples. Signal and background estimates for the dilepton 
channel are summarized in Table IIIII Uncertainties are 
due to integrated luminosity, MC sample statistics and 
fake rates. 



G. MC simulation 

The signal (tt) MC simulation is modeled by PYTHIA 
version 6.216 [111, with herwig version 6.510 29] used 
as a cross-check. Most background kinematics are esti- 
mated from MC samples. The diboson backgrounds are 
modeled with pythia version 6.216 and the W + jets 
and Drell-Yan+jets backgrounds are modeled by ALP- 
GEN version 2.10' [s^l, with jet fragmentation modeled 
by PYTHIA version 6.325 A matching scheme [3l| is 
used to ensure that there is no double-counting of phase 
space in background events, as it is otherwise possible 
for events with hard hadronic shower evolution to give 
states already described by events at the matrix element 
level. The W + jets and Drell-Yan+jets background MC 
is divided into 7i-parton samples, where n refers to the 
total number of partons (quarks or gluons), including 
heavy flavor. The samples are combined according to the 
cross sections reported by ALPGEN, accounting for pos- 
sibly different efficiencies for the samples to pass event 
selection. Similarly, events with heavy flavor after frag- 
mentation are checked to ensure no double-counting of 
phase space across samples with different flavor types at 
the matrix element level. Double-counting can occur if 
samples generated at the matrix element level with light- 
flavor partons produce charm or bottom quark pairs in 
the parton shower. We remove such events unless both 
heavy flavor partons are within the same jet. In addition 
the heavy flavor quark pairs generated at the matrix el- 
ement level can enter the same jet effectively reducing 
heavy jet multiplicity, therefore we also remove events 
of this type. Depending on the multiplicity of light and 
heavy flavor final state partons we remove a fraction of 
events between a few percent and approximately 20%. 
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Electroweak production of single top quarks in both 
the s- and ^-channels contributes very few events to our 
sample. These events are treated as background, and 
are modeled using a fixed mass of Mtop = 175.0 GeV/c^. 
Single top quark events are generated by madevent 33] ; 
fragmentation is modeled with pythia version 6.409 33] ■ 

To model multiple proton- antiproton interactions oc- 
curring in a single bunch crossing, we add interactions 
where no partons with high transverse momenta are pro- 
duced to the events simulated for each process. Those 
minimum bias collisions are simulated with pythia ver- 
sion 6.216 . The number of minimum bias interactions 
added to a given event is equal to the expected number 
of pp interactions, which depends on the instantaneous 
luminosity profile of the data run of the event. The in- 
stantaneous luminosity profile is matched between MC 
samples and data only for the first 1.2 fb^^ of integrated 
luminosity. This incorrect model is a source of bias of 
0.4 GeV/c^ for the dilepton-only fit. No bias is present 
in the fitted top quark mass in the combined and lepton 
-I- jets only fits, however a bias of 0.04 CTc is present in the 
fitted AjES in both measurements. The bias in the dilep- 
ton measurement is higher due to the particular choice 
of observables used to make the measurement and also 
due to the fact that the in situ calibration absorbs the 
bias on mass and converts it into bias on A jes ■ These bi- 
ases were found in studies of MC samples with increased 
number of pp interactions and are corrected for in quoted 
results. 



IV. EVENT RECONSTRUCTION 

After selection and reconstruction of event parameters 
of physics interest, the data are processed to form esti- 
mators for the top quark mass. Simply forming invariant 
masses is not possible, as there arc two top quarks per 
event. Events in the lepton-|-jets channel have many pos- 
sible assignments of jets to the quarks, each of which give 
different reconstructed top quark masses. Each dilepton 
event has two undetected neutrinos, resulting in under- 
constrained kinematics. Both channels must account for 
the possibility that, due to radiation effects, jets in the 
detector may not correspond to quarks from the hard 
scattering. Both topologies also contain non-negligible 
backgrounds. We approach these problems by construct- 
ing quantities strongly correlated to the top quark mass 
and comparing the data to MC predictions that include 
all of the above effects. For the lepton-|-jets channel, we 
additionally account for the unique and dense environ- 
ment of tt events by applying jet corrections specific to 
the ti events. 

A set of generated distributions for a particular top 
quark mass is referred to as a template. A template is 
then a probability density function for a set of observ- 
ables. Our measurement of the top quark mass is then a 
determination of the most likely parent template for the 
data. Further complicating the analysis, however, is the 



strong correlation between the Ajes in the detector and 
quantities sensitive to the top quark mass, including the 
top quark mass estimators. As explained in Sec. IIIIBl 
scaling measured jet energies back to original parton en- 
ergies is a difficult task, and any uncertainty on the JES 
directly translates to a systematic uncertainty on Mtop- 
To reduce this effect, we introduce a second template in 
the lepton-t-jets channel that uses the hadronic decay of 
the W boson to make an in situ measurement of Ajes 
that can be applied to all jets in the event sample, includ- 
ing those from b quarks and those in the dilepton channel. 
The narrow width of the W makes its dijet mass {rrijj) 
a good estimator for Ajes. 



A. Top-quark specific corrections 



The jet corrections described in Sec. IIII Bl are generic 
algorithms derived for application in all high-energy CDF 
analyses. As such, they miss out on several key features 
of tt events in the lepton-|-jets channel. The generic jet 
corrections assume fiat px spectra for all jets. The bias 
for the pt spectra expected from the physics process 
under consideration may be corrected specifically. Top- 
antitop events have two different, non-flat px spectra for 
the W decay jets and the b jets. The generic jet cor- 
rections also do not account for differences between jets 
coming from b quarks and jets coming from light-ffavor 
quarks. To account for all these effects, we derive jet cor- 
rections specific to the lepton-|-jets tt environment from 
MC simulations. After event selection, jets arc corrected 
to the particle jet level. We separate the corrected jets 
by fiavor (whether they came from a b quark or a light 
quark), and then into different j] and px bins. The top- 
specific corrections are derived such that the px of a jet 
corresponds to the most probable value of the quark pro- 
ducing the jet. 

The top-specific corrections can be over 50% for low-px 
jets in the central region, and slightly negative for high- 
px jets. We apply the corrections to the px of the jet. 
We assume the direction of the jet to be well measured, 
and do not apply correction to the jet angles {ri and (/>). 
To obtain the jet energy, the mass of jets assumed to 
come from b quarks is fixed to 5.0 GeV/c^ and the mass 
of jets assumed to come from light quarks is fixed to 
0.5 GeV/c^, though the mass effects are small compared 
to typical jet energies in tt events. The top-specific cor- 
rections also provide the resolution on jet energy, once 
again separately for the two fiavors and as a function of 
Px and 77. The resolution is worse than 20% for low-px 
jets, and better than 10% for high-pT jets. 



B. Lepton+jets reconstruction 

The lepton-|-jets decay channel gives overconstrained 
kinematics for the tt system. Detailed information on 
the lepton-fjets kinematic fitter can be found in Ref. [3]. 
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The minimization package MINUIT [3J| is used to mini- 
mize a x^-like function for the overconstrained kinematic 
system: 



i—i,4 jets 



i,fit i,mcas\2 



(Pi - Pi 



E 



^ w 



^ w 



(3) 



The first term constrains the pt of the lepton and the 
4 jets in the event to their measured values, within their 
uncertainties ct^. The unclustered energy ([/) is the en- 
ergy in the calorimeter not associated with the primary 
lepton or one of the four leading jets. The second term 
constrains the x and y components of the unclustered 
energy {Ux = ?7sin6'cos0, Uy = Usm9sm(j>) in the de- 
tector close to their measured values within uncertainties 
aj. The third term in the expression constrains the 
dijet mass of the two jets assigned as W decay daughters 
to the well measured W mass within the W boson decay 
width. The fourth term similarly constrains the invari- 
ant mass of the leptonic W decay daughters. The last 
two terms constrain the invariant masses of the three- 
body top decay daughters to be consistent within the 
top quark decay width of 1.5 GeV. The value of ml°'^° 
is a free parameter in the fit, and is taken as the recon- 
structed mass used in the templates. 

The neutrino transverse momentum is not a direct pa- 
rameter in the minimization, but is instead related to 
the unclustered energy, and is calculated at every stage 
of the minimization process: 



Px.yH 




(lepton) -I- Ux 



(4) 



The longitudinal component of the neutrino momentum 
is a free parameter that is effectively determined by the 
constraint on the invariant mass of the leptonic W . 

With the assumption that the leading (most energetic) 
four jets in the detector come from the four final quarks 
at the hard scatter level, there are 12 possible assign- 
ments of jets to quarks. The minimization is performed 
for each assignment, with ml'^'^° taken from the assign- 
ment that yields the lowest x^- Events with the lowest 
> 9.0 are removed from the sample to reject poorly 
reconstructed events not fitting the tt hypothesis. The 
cut was optimized for expected statistical precision how- 
ever we find that there is no strong dependence of the 
expected precision on the value of the cut. Identifying b 
jets reduces the number of combinations since tagged jets 
are assigned only to final state b quarks. In rare events 
with more than two tags among the leading four jets, 
only two tags of highest Et are assigned to the b quarks 
and the additional tags are ignored. 



The calculation of the dijet mass rrijj is independent of 
the above minimization procedure to derive mj"'^". Given 
a pair of jets, a simple invariant mass is calculated from 
the jet four-vectors; in particular, the W mass constraint 
of the kinematic fitter is not applied. There are multiple 
ways to choose two jets among the four or more jets in tt 
events. Tagged jets are assumed to come from final state 
6 quarks. Additionally, the two jets from the hadronic W 
decay daughters are assumed to be among the leading 
four jets. For two-tag events, there is only one choice 
for the jet pair to be associated to the W boson. For 
1-tag events, there are three possible dijet masses to be 
made from the 3 non-tagged leading jets; we pick the 
single dijet mass closest to the well known W mass. This 
sculpts the distribution, but is the choice most likely to 
be correct in selecting the two jets from the W decay 
daughters, and was found to give the best sensitivity to 
AjEs- We correct the jets using the light quark top- 
specific corrections. 



C. Lepton+jets template results 

We process MC samples with different values of Mtop 
with full detector simulation [35] and event selection. 
The kinematic fitter is applied to each event, giving the 
ml'^'^° templates shown in Fig. [5] (a) and (c). Though 
the peak of the templates depends strongly on Mtop, the 
reconstruction is not perfect, and ml°'^° only gives an es- 
timate for Mtop. The large tails in the templates are a 
result of incorrect jet-quark assignments. The 2-tag sub- 
sample, with fewer jet-parton assignments, has narrower 
templates, and therefore has more sensitivity to the top 
quark mass. Templates for rrijj masses for three different 
values of A jes are shown in Fig. [2] (b) and (d) . The 1- 
tag rrijj templates are narrower than the 2-tag templates 
due to sculpting of the distributions. The sculpting also 
yields smaller shifts in the 1-tag template then in the 
2-tag template, as the Ajes varies. 



D. Dilepton templates 

In the dilepton channel, the measured quantities and 
assumptions on the masses of particles in the decay cas- 
cade do not provide enough constraints to reconstruct 
the four-vectors of the top quarks. Instead, we form a 
reconstructed top quark mass (mf^^) for each dilep- 
ton candidate using the Neutrino Weighting Algorithm 
[1, Q. The algorithm assigns a weight to the event as 
a function of top quark mass. A top quark mass scan 
is performed in the range 80-380 GeV/c^, and the value 
yielding the maximum weight is selected as mf^^. The 
two most energetic jets in an event are considered to have 
originated from the b quarks, giving two possible jet- 
quark assignments. For each assumed top quark mass 
and jet-quark assignment, we integrate numerically over 
the possible pseudorapidities of the two neutrinos. The 
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FIG. 2: Template distributions for MC events passing the lepton + jets selection. Shown are the 1-tag m\^'^° (a) and mjj (b) 
distributions, and the 2-tag m\'"^° (c) and nijj (d) distributions. The m\'''^° distributions are plotted using events with three 
values of Mtop and with the nominal Ajes = 0.0 ctc. The rrijj distributions are plotted using events with three values of Ajes 
and with Mtop = 170 GeV/c^ 



distribution of the neutrino pseudorapidity is assumed to 
be Gaussian around zero, with a width of 1.0 obtained 
from PYTHIA. Given a neutrino pseudorapidity, we can 
solve for its transverse momentum. Up to two solutions 
are possible for each of the neutrino and antineutrino 
transverse momenta. For each of the four solutions, we 
compare the total momentum in the x and y directions 
carried by the neutrinos to the measured x and y com- 
ponents of the We calculate a weight that is the 
product of two Gaussians, one each for the x and y direc- 
tions, of the difference between the measured 6t in that 
direction and the sum of the momentum components of 
the two neutrinos. We use a Gaussian width of 19 GeV, 
which is optimized using MC tt events. The four weights 
are added to form the integrand. Note that we do not 
account for resolution effects in measurements of the jets 
and leptons. The two integrals corresponding to the two 
jet-6 quark assignments are added to form the top-mass- 
dependent weight. In the calculation of the transverse 
momenta of the neutrinos, jets are corrected using the 



generic jet corrections; applying the top-specific correc- 
tions of Sec. lIV Al was not found to improve the resolution 
on the reconstructed top quark mass in this channel. Fig- 
ure [3] (a) and (c) shows the output of the algorithm from 
fully simulated MC events with different input masses. 

The momenta of the decay products of the tt pair are 
directly correlated to the invariant masses of the top 
quarks. We therefore use another variable, the Ht [cf. 
Section IIII F[) , to improve the precision of the measure- 
ment, quark mass. Figure[3](b) and (d) shows Ht distri- 
butions for different top quark masses. The correlation 
coefficient between m^^^ and Ht is '-^ 40% in signal 
events and ~ 60% in background events. 



V. MASS FITTING 

The distributions of the observables defined in Sec. IIVI 
are used to determine simultaneously the two param- 
eters Mtop and Ajes- For this, we need to know 
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(c) Tagged rrit (d) Tagged Ht 

FIG. 3: Template distributions for MC events passing the dilepton selection. Shown are the 0-tag mf^^ (a) and Ht (b) 
distributions, and the tagged mf^^ (c) and Ht (d) distributions. The distributions are plotted using events with three values 
of Mtop and with the nominal Ajes ~ 0.0 a^. 



P(x, y; Mtop, AjEs): the probability of observing a par- 
ticular pair of values of the mass-sensitive parameters 
{x,y), given some Mtop and Ajes- The observables x 
and y correspond to m™'^° and rrijj for lepton -I- jets 
events and to m^^^ and Ht for dilepton events. When 
Mtop and Ajes are fixed, the resulting P(x, y) should 
be a normalized probability density function (PDF) over 
the two-dimensional space of the observables. The PDFs 
must be determined separately for signal and background 
events in each subsample (e.g. 1-tag lepton -I- jets events). 
The background probabilities do not depend on Mtop. 

Inaccuracies in these families of PDFs lead to biases in 
the final measurement that can be difficult to uncover or 
to characterize. Therefore, in order to achieve a precision 
measurement of Mtop , it is essential to make a robust de- 
termination of P for each class of events. We accomplish 
this in two steps: First, at discrete values of Mtop and 
Ajes I we estimate the two-dimensional PDFs for the ob- 
servables from large samples of MC events using kernel 
density estimation (KDE), described in Sec. IV Al Then 
we smooth and interpolate to find PDFs for arbitrary val- 



ues of Mtop and Ajes using local polynomial smoothing 
(LPS), described in Sec.lTBl 

The resulting probabilities are used in a combined like- 
lihood fit (Sec. IV Cp to measure the top quark mass. We 
run a rigorous set of checks to validate the analysis ma- 
chinery and calibrate the final result using events from 
MC simulation; these checks are described in Sec. IV Dl 

A. Kernel Density Estimation 

Previous template-based measurements of the top 
quark mass j36i | used arbitrary functional forms to fit 
parameterized PDFs from histograms of the observables. 
It is difficult to extend such parameterizations to two 
dimensions in observables while properly accounting for 
correlations between m\'^'^° and mjj, or m^^^ and Ht- 
In the lepton -I- jets channel, these correlations can lead 
to a bias of several hundred MeV, if not properly taken 
into account; in the dilepton channel, the correlations 
are larger and can make such a measurement with two 
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observables impossible. Using a functional form with a 
large number of parameters can also result in fits that are 
unstable with respect to small changes in the histograms 
or even in the parameter initialization. 

This measurement takes a different approach based on 
KDE to form PDFs in two observables without any as- 
sumption about the functional form. Some useful intro- 
ductions to KDE can be found in Refs. [32, i38i, i39|. In 
KDE, the probability for an event with observable x is 
given by a linear sum of contributions from all entries in 
the MC sample. For a one-dimensional distribution, this 
probability is given by: 



n 

nh ^ 



(5) 



where f{x) is the probability to observe x given, as an 
example, a tt MC sample with known Mtop and Ajes- 
The sample has n entries, with values of the observable 
given by Xi. The kernel function K is a normalized func- 
tion that adds less probability to a measurement at x as 
its distance from Xi increases. The smoothing parame- 
ter h (sometimes called the bandwidth) is a number that 
determines the width of the kernel. Larger values of h 
smooth out the contribution to the kernel density esti- 
mate and give more weight at x farther from Xi. Smaller 
values of h provide less bias to the kernel density esti- 
mate, but are more sensitive to statistical fluctuations. 
We use an Epanechnikov kernel, defined as: 



for 1^1 < 1, 
otherwise. 



(6) 



so that only events with |a; — x^j < h contribute to f{x). 

We use an adaptive KDE method in which the value of 
h is replaced by hi so that the amount of smoothing ap- 
plied to the ith event depends on the value of f{xi) [40| . 
We run a first pass of kernel density estimation with con- 
stant h. This pilot kernel density estimate is then used 
in a second round of KDE to determine the individual 
hi, with hi cx f{xi)^^-^. In the peak of the distribu- 
tions, where there are more events, we use small values 
of hi to capture as much shape information as possible. 
In the tails of the distribution, where there are fewer 
events and the kernel density estimates are sensitive to 
statistical fluctuations, a larger value of hi is used. The 
overall scale of h is set by the root mean square (RMS) of 
the distribution and by the number of entries in the MC 
sample; larger (smaller) smoothin g is used when fewer 
(more) events are available [33,[3l,l4lj. If the smoothing 
parameters get too large in the tails of the distribution, 
the kernel density estimates can become non-local, and 
a point at Xi can contribute weight to an estimate at a 
distant x. Following Ref. [l^l, we guard against this by 
not allowing hi to get too large: 



(7) 



where /iq is the minimum adaptive bandwidth, which oc- 
curs in the peak of the distribution. 

KDE is extended to two dimensions by multiplying two 
kernels together [H, : 



n 

-E 



i=l 



K 



X K 



hi = min(/ii, v 10 • /iq) 



(8) 

Note that the smoothing parameters for the two variables 
do not have to be identical. Typical values of h for ker- 
nel density estimates in the signal are 10-12 GeV/c^ for 
m\'"'°, 3-6 GeV/c^ for m^j, 15-20 GeV/c^ for m^^^ and 
45-55 GeV for Ht- For background kernel density esti- 
mates, these numbers are slightly larger, as the number 
of events passing all the cuts is smaller and the templates 
are wider. 

The kernels in Eq. ^ know nothing about the bound- 
aries of the templates. Mathematically, the density func- 
tions can take on any real numbers given large enough 
smoothing, even though kinematic requirements and en- 
ergy conservation limit possible values of the observables. 
When the probability density extends beyond such a 
limit, where the data are not found, the normalization 
condition of the kernels does not hold. To enforce unit 
normalization, we explicitly force hard boundaries and 
reject events in the tails of the distribution both from the 
MC templates and the data, typically removing 1-2% of 
signal events and a slightly larger fraction of background 
events. When kernel density estimates are calculated, 
we check that each of the individual kernels is normal- 
ized within the boundaries. If a kernel is not normalized 
and leaks probability outside the boundaries, it is renor- 
malized such that it contributes unit weight inside the 
boundaries. 

Using f{x,y) from ([5]) as P{x,y), we can scan val- 
ues of X and y to visualize the two-dimensional PDF for 
fixed Mtop and Ajes ■ Figure 3] shows two-dimensional 
kernel density estimates for lepton -I- jets and dilepton 
signal events, given Mtop = 170 GeV/c^ and Ajes = 0.0. 
Figure [5] shows the estimates for background events at 
Ajes = 0.0. The background kernel density estimates 
are derived separately for the individual contributions to 
the background model, taking into account the sample 
sizes and RMS values, and are then combined with the 
appropriate weights. Note that since jet energy scale 
shifts would affect all data events in a similar way, data- 
derived background templates do not depend on Ajes. 
The dilepton tagged background contains multiple peaks 
that come from the fake background. The data used to 
model the fakes, which comprise half of the total tagged 
background in the dilepton channel, contain very few 
events, and also peak at different locations than the other 
backgrounds. 
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FIG. 4: Kernel density estimates at Mtop ~ ITO GeV/c^ and Ajbs = 0.0 for lepton + jets 1-tag (a) and 2-tag (b) events, and 
for dilepton untagged (c) and tagged (d) events. 



B. Local Polynomial Smoothing 

The PDF families P(a;, y; Mtop, Ajes) are defined on 
the continuous parameters Mtop and Ajes- But the sig- 
nal MC samples are produced at discrete values of Mtop 
and Ajes, and the background MC samples are pro- 
duced at discrete values of Ajes, so KDE is capable of 
producing PDFs only at discrete points of the param- 
eter space. In addition, statistical fluctuations in the 
kernel density estimates are correlated for events with 
similar observables, so it is useful to smooth out the 
PDFs before the likehhood fit. To obtain PDFs that are 
smoothly and continuously varying as a function of Mtop 
and Ajes without assuming Gaussian likelihoods, we em- 
ploy a technique known as local polynomial smoothing 
(LPS) nil, described briefly below. 

LPS locally approximates the value of the PDF with 
a second-order polynomial. The expansion uses the esti- 
mates from KDE, but gives more weight to MC samples 
in a nearby region of Ajes (and Mtop , if we are smooth- 



ing out the signal probabilities) . We look for an estimate 
P(a) for the true value of the function P{(y.), where we 
have omitted the values of the observables {x,y) from 
the arguments of the function. The quantity a is a two- 
dimensional vector (Mtop, Ajes) in the case of the sig- 
nal probability function or the scalar Ajes for the back- 
ground. Kernel density estimation provides estimates Yfc 
for the values of P{a.k) at a number of points a.^. We 
assume that Yfc are unbiased estimators of the true prob- 
ability values P{ak) with the same variance. 

A second order expansion of the function P for points 
t in the neighborhood of a can be written as: 

P{t) = {c,Fit-a)) (9) 

where the angle bracket denotes an inner product. The 
coefficients of the expansion are given by the compo- 
nents of the vector c. The quantity is a vector of 
basis functions for second order polynomials. For a two- 
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FIG. 5: Kernel density estimates at Ajes = 0.0 for lepton + jets 1-tag (a) and 2-tag (b) background events, and for dilepton 
untagged (c) and tagged (d) events. 



dimensional v, F{v) is: 



Fiv) 



( 1 \ 



(10) 



If ?; is a scalar, Fiv) reduces to a 3-component vector. 

To evaluate P{ct), we minimize the criterion given by 
Eq. (fTTj) with respect to c. In other words, we find the 
second order expansion of P around a that best matches 
the estimates Yk at points 0;^: 



(11) 



where the weight of each estimate is given by the factor 



w{otk) = W{uk), with: 



Uk 
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E 



OLd - OLk,d 



(12) 



The sum in Eq. (|12ll runs over the components of the 
vectors, and is the dimensionality of the parameter 
space. We use W = W{u) = (1 — \u\^)^ for |m| < 1 and 
W{u) = otherwise. This gives a smoothly decreasing 
weight to the estimates Yk obtained at points far away 
from the evaluation point a. The constants control 
the amount of smoothing, larger values of give more 
weight to the estimates farther away from the point a. 

For this analysis, signal MC samples are generated 
at 76 mass points with Mtop ranging from 120 to 
240 GeV/c^. The spacing between mass points is small 
(0.5 GeV/c^) in the region of interest (165-185 GeV/c^), 
and gets larger in the tails of the grid. Each signal 
MC sample and MC-based background is processed us- 
ing 29 different values of Ajes from —3.0 (Tc to -1-3.0 ctc. 
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Near the range of interest of nominal Ajes, the spac- 
ing is 0.2 (Tc- We set the constants hd of Eq. based 
on the performance of the analysis in terms of expected 
precision and biases. Wc choose /iMtop = 10.0 GeV/c^ 
for signal PDF smoothing in lepton + jets events and 
/iMtop — 15.0 GeV/c^ in the dilepton events. For both 
categories, /iajes = 0.8 CTc- The background MC sam- 
ples have smaller statistics and increased jitter, as many 
events are selected near jet energy thresholds and move 
in and out of the sample as Ajgg varies. To compensate 
for these effects, we choose the larger /iajes = 3.0 cTc for 
background smoothing. 



C. Likelihood fit 

We compare the two dimensional distributions of the 
observables in the data with the signal and background 
PDFs in an unbinned extended maximum likelihood 
fit [i^. The most important parameters of the fit are 
the mass of the top quark (Mtop) and the deviation from 
the nominal jet energy scale (Ajes)- Each subsample 
gives two additional parameters: the expected number 
of signal events (jig) and the expected number of back- 
ground events (rib) in the subsample. The likelihood form 
is given by: 

C — e 2 

cptoii+jcts,l-tag 
^ •^lcpton+jcts,2-tag 



^ ^dilepton, non-tagged 
^ -^dilepton, tagged 



(13) 



where the first term in the product constrains the mea- 
surement of the JES to its nominal value of within the 
uncertainty of 1 CTc- Each of the subsequent terms corre- 
sponds to one subsample, and is given by: 



sample 



shape 



X L 



bg 



(14) 



The term in the likelihood most critical to the mass 
measurement is the extended maximum likelihood shape 
term: 



shape 



N 
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n 



m 

nsPs{xi,yi; Mtop, Ajes) + nbPb{xi,yi] Ajes) 



Us + Jib 



(15) 



where the product runs over all events in a given 
subsample. The observables {xi,yi) in the ith event 
are {m\^'^° ^rajj) in the lepton -I- jets channel and 



(m 



NWA 



Ht) in the dilepton channel. The quantities Ps 



and Pb designate the signal and background PDFs as de- 
termined by KDE and LPS. To improve the precision of 



the measurement, we apply a Gaussian constraint to the 
expected number of background events: 



(16) 



where Ub^ is the a priori estimate for the expected num- 
ber of background events, and f7„^^ is the uncertainty on 
the estimate. Both sets of numbers are given in Table HIl 
and Table UTTl 

We minimize the negative logarithm of the likelihood 
with respect to all 10 parameters using minuit. The un- 
certainty on Mtop and Ajes is found by searching for 
the points where the negative logarithm of the likelihood 
minimized with respect to all other parameters deviates 
by i from the minimum. The uncertainty on the top 
quark measurement obtained this way includes the sta- 
tistical uncertainty as well as the systematic uncertainty 
due to allowed variations in the jet energy scale and the 
background estimates. We scale the uncertainty on the 
top quark mass by the pull width as obtained in Sec. lVDl 
For the dilepton-only cross-check, we fix the value of the 
Ajes parameter to Uc and perform a fit only for Mtop, 
as the dilepton channel has no power to resolve Mtop and 
Ajes simultaneously. 



D. Metiiod checlc 

We test the likelihood procedure using large numbers 
of MC simulated experiments, each of which is generated 
for a specific value of Mtop and Ajes- In each experi- 
ment we select the number of background events from a 
Poisson distribution with a mean equal to the expected 
number of background events in the sample. The num- 
ber of signal events is selected from a Poisson distribu- 
tion with a mean equal to the expected number of sig- 
nal events assuming a tt pair production cross section 
of 6.7 pb. The signal events are drawn at random from 
a MC sample generated at a given Mtop and Ajes- The 
background events are drawn from the entire background 
sample with probabilities corresponding to the different 
background sources and the weights of individual events 
given by the model. Individual event rates can vary based 
for example on mistag probability of jets in a given event. 
An event (signal or background) can be drawn arbitrary 
number of times from a given MC sample so that differ- 
ent MC experiments can share simulated data. We refer 
to this method of drawing events from MC samples as 
"drawing with replacement" . Once the simulated data 
is constructed we perform a maximum likelihood fit as 
described in the previous sections. We fluctuate the con- 
straint on JES and the constraint on the expected num- 
ber of background events in each subsample to reflect the 
possible deviation of those parameters from their a pri- 
ori expected values. The fluctuations on the constraints 
are applied to estimate the effect of our limited knowl- 
edge about the nuisance parameters on the top quark 
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mass measurement. The jet energy scale constraint in 
Eq. (fT51) is replaced by: 

e (17) 

where Ajes^ is randomly selected in each MC simulated 
experiment from a Gaussian with mean corresponding 
to the AjEs being tested and unit width. Similarly, the 
quantity in Eq. (|16p is replaced in each MC simulated 
experiment by a value drawn from a Gaussian with mean 
of nfcg and width cr„j_^ . 

The likelihood fit should, on average, return the value 
of the top quark mass used to generate the MC simulated 
experiments. Figure [H] shows the average residual (de- 
viation from expectation for an unbiased measurement) 
from 3000 MC simulated experiments for a range of Mtop- 
The fit to a constant shows no bias for the combined 
and lepton + jets-only fits, and a small positive bias for 
the dilepton-only fit. This bias does not warrant a cor- 
rection since it is small in comparison to the expected 
uncertainty, and has a probability of 9% for a purely sta- 
tistical fluctuation to generate the observed shift. 

The statistical uncertainty on the measurement is ex- 
tracted from the data. To test whether the error estimate 
is sound, for each MC simulated experiment we calculate 
the pull, defined as a ratio of the residual to the un- 
certainty reported by MINUIT. If the residual is positive 
(negative) we use the negative (positive) error in the ra- 
tio. The width of the pull distribution for Mtop is shown 
in Fig. [T] The average pull width is larger than 1.0 for 
the combined and lepton -I- jets-only measurements due 
to the finite number of events in the two-dimensional fits. 
The pull width correction is 3% for these measurements, 
and thus we increase the reported uncertainty in the data 
by this amount. 

The residual and pull width for the Ajes parameter 
are also investigated using the MC simulated experiment 
ensembles. Both the combined fit and the lepton -I- jets- 
only fit show a negative bias of < 0.02 CTc (with marginal 
statistical significance) and a pull width of ^ 1.04. Since 
Ajes is a nuisance parameter, whose precise value is not 
as important as its effect on the top quark mass mea- 
surement, we do not correct the Ajes value measured in 
data for this bias. 

As noted above, when ensembles of simulated data are 
constructed the events are drawn with replacement. If 
we were to draw the events from MC samples without 
replacement such that no two MC simulated experiments 
were to share events, we would have only ~ 100 MC sim- 
ulated experiments for each ti sample. Drawing events 
with replacement allows us to perform an arbitrarily large 
number of MC simulated experiments, fully exploring the 
possible combinations of events in order to check our ma- 
chinery for possible biases. To evaluate the uncertainties 
on statistics such as the residuals and pull widths, we 
employ the bootstrap technique [i^, 113] . In each boot- 
strap ensemble, we draw events from the signal MC sam- 
ple with replacement until we reach the same number of 



events as in the original sample. We then run 3000 MC 
simulated experiments using this bootstrapped sample in 
place of the original sample. We repeat the above proce- 
dure 60 times. For each of the bootstrap ensembles, we 
calculate the desired statistic. The RMS of the statis- 
tic in question from the 60 bootstrap ensembles is taken 
as the uncertainty on the statistic. For example the un- 
certainty on the residual of fitted top quark mass in a 
typical MC sample is 0.4 GeV/c^ for the dilepton fit and 
0.2 GeV for the lepton -I- jets and combined fits 



VI. RESULTS ON DATA 

The likelihood fit when applied to the data yields 
Mtop = 171.9 ± 1.7 (stat. + JES) GeV/c^. The 
lepton -I- jets-only fit yields Mtop = 171.8 ± 
1.9 (stat. -I- JES) GeV/c2. The dilepton-only fit, 
which does not include an in situ Ajes measure- 
ment but instead fixes Ajes to 0.0 CTc, yields Mtop — 
171.2 (stat.) GeV/c^. The combined fit returns 

Ajes = -0.17 ± 0.35 (stat. + Mtop) ctc, and the lepton 
+ jets fit returns Ajes = -0.12±0.37 (stat. + Mtop) ^c. 
The results above have been corrected for the pull width 
and high instantaneous luminosity effects. Results from 
the combined fit, including fitted numbers of signal and 
background events for each subsample, are summarized 
in Table [TVl The dilepton-only fit and lepton -f jets-only 
fit both return Mtop value lower than the Mtop measured 
in both channels simultaneously. This is due to the in situ 
JES calibration extracted from the lepton -I- jets channel 
events being applied to the dilepton channel data in the 
combined fit. 

The log-likelihood contours for the combined measure- 
ment are shown in Fig. [51 The one-dimensional log- 
likelihood for the dilepton-only measurement and the 
log-likelihood contours for the lepton + jets-only mea- 
surement are shown in Fig. [H] and Fig. 1101 respectively. 
Figure [TT] shows the one-dimensional lepton -I- jets data 
with the best-fit one-dimensional signal and background 
distributions overlaid on top. Figure [T^ shows the distri- 
butions for the dilepton data. Using the observed number 
of events in data and the background expectations, 10% 
of MC experiments have a smaller error than the value 
measured in the combined fit. The p- value for the lepton 
+ jets-only fit is 21%; the value for the dilepton-only fit 
is 14%. 

We fit without the a priori JES and background con- 
straints and measure the same Mtop = 171.9±1.7 (stat.-l- 
JES) GeV/c^, showing that these priors do not signif- 
icantly affect our result. We also fit separately in the 
several individual subsamples: the first 1 fb~^ and last 
0.9 fb^^ of data, electron and muon events in the lepton 
-I- jets-only fit, and different lepton pair type events in 
the dilepton-only fit. In addition we quote the top quark 
mass fitted in lepton -I- jets and dilepton subsamples sep- 
arated by b tagging multiplicity. The results are consis- 
tent across cross-checks and are summarized in Table [V] 
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FIG. 6: Cliecks for bias in the fitted top quark mass for a) lepton + jets only fit, b) dilepton only fit and c) combined fit. 



TABLE IV: The input constraints and fitted values are given 
for all free parameters in the combined likelihood fit. LJ refers 
to lepton + jets subsamples and DIL refers to dilepton sub- 
samples. Con. refers to the constraint used in the likelihood. 



Catej 


;ory 


LJ 2-tag 


LJ 1-tag DIL tagged 


DIL 0-tag 


Mtop 


con. 
fit 




None 

171.9 ± 1.7 GeV/c^ 




JES 


con. 
fit 




0.0 ± 1.0 ac 
-0.17 ±0.35 (Jc 




Us 


con. 




None 






fit 


96.4 


184.1 till 56.9 t.7,1 


43.8 11°,' 


rib 


con. 


4.2 ± 1.9 


42.7 ±12.5 3.9 ±1.0 


44.3 ± 7.0 




fit 


3.4 ± 1.9 


47.4 3.9 ±1.0 


41.5 til 



The results are not corrected for pull width effects and 
bias due to incorrect instantaneous luminosity profile of 
the MC samples. With the exception of the fit without 
the a priori JES constraint, all cross-checks include the 
JES prior. The fits in the four subsamples separated by b 
tagging multiplicity include the a priori background con- 
straints; all other cross-checks do not include background 
constraints. 



VII. SYSTEMATIC UNCERTAINTIES 

We examine a variety of effects that could affect our 
measurement by comparing MC simulated experiments 
in which we change systematic parameters within their 
uncertainties. As a single nuisance parameter, the mea- 
sured A JES does not fully capture the complexities of jet 
energy scale uncertainties, particularly those with differ- 
ent T] and pt dependence. Fitting for the global JES 
removes most of these effects, but not all of them. In 
order to estimate the total residual JES uncertainty, we 
vary JES parameters within their uncertainties in both 
signal and background MC generated data and measure 
resulting shifts in Mtop- We also conduct MC simulated 
experiments where we assume JES uncertainties are not 
fully correlated between jets of different momenta. So 
as not to bias the results, we remove the JES prior for 
these experiments. For the dilepton-only measurement, 
which has no in situ calibration, these systematics dom- 
inate. To form a b jet energy scale systematic wc replace 
the default parameters of the Bowler fragmentation func- 
tion [48| in the PYTHIA simulation with the parameters 
obtained by the D0 collaboration in a PYTHIA tune to the 
LEP and SLD data [i^. We also vary the semileptonic 
branching fractions of b and c quarks within their uncer- 
tainties given by [s^l and [l|. The calorimeter response 
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FIG. 7: Width of the pull distribution for the fitted Mtop for the a) lepton + jets only fit, b) dilepton only fit and c) combined 
fit. 
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FIG. 8: Negative log-likelihood contours for the combined fit. 
The minimum is indicated by the 'x' and corresponds to the 
most probable top quark mass and Ajes, given the data. The 
contours are drawn at values of 0.5, 2.0 and 4.5 of the increase 
of the log-likelihood from the minimum value. These curves 
correspond to the 1,2 and 3 a uncertainty on the measurement 
of the top quark mass. 



FIG. 9; One-dimensional log-likelihood for the dilepton-only 
fit. 



to b jets is varied to capture differences in absolute jet 
energy scale uncertainties for light flavor and b quarks. 
Effects due to uncertain modeling of initial-state gluon 
radiation (ISR) and final-state gluon radiation (FSR) are 
studied by extrapolating uncertainties in the px of Drell- 
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FIG. 10: Negative log-likelihood contours for the lepton + 
jets-only fit. The minimum is indicated by the 'x' and corre- 
sponds to the most probable top quark mass and Ajes, given 
the data. The contours are drawn at values of 0.5, 2.0 and 4.5 
of the increase of the log-likelihood from the minimum value. 
These curves correspond to the 1,2 and 3 a uncertainty on 
the measurement of the top quark mass. 



TABLE V; Cross-checks on the data. LJ refers to the lepton 
-I- jets-only fit, DIL refers to the dilepton-only fit, and Combo 
refers to the combined fit. All numbers are uncorrected for 
pull width effects and bias due to incorrect instantaneous lu- 
minosity profile of the MC samples. For the dilepton-only fits, 
Ajes is fixed to 0.0 CTc. 



Fit type 


Sample 


Mtop (GeV/c^) 


Ajes (o-c) 




Combo 


171.9 ± 1.7 


-0.12 ±0.34 


Nominal 


LJ 
DIL 


171.8 ± 1.8 
171.6 t^ j 


-0.09 ±0.36 




Combo 


171.9 ± 1.7 


"J-J^^ -0.37 


No JES prior 


LJ 


171.8 ± 1.9 


-0 11 +"■3^ 




DIL 


171.6 tli 




Combo 


171.9 ± 1.7 


— 1 1 +"'^^ 


No bkgd prior 


LJ 
DIL 


171.8 ± 1.8 
171.5 ±3.4 


-0.06 ±0.36 


1-tag LJ 


LJ 


169.1 lis 


-0 17 

-0.57 


2-tag LJ 


LJ 


173.6 lal 


20 +" *'^ 

"•^U -0.50 


0-tag DIL 


DIL 


170.1 1?;^ 




Tagged DIL 


DIL 


172.2 




e only LJ 


LJ 


172.2 ±2.7 


-0.09 ±0.51 


/i only LJ 


LJ 


171.3 t'tt 


-0 04 

^■'-'^ -0.47 


ee only DIL 


DIL 


169.0 ±8.0 




e/i only DIL 


DIL 


173.6 tl-'i 




Hfi only DIL 


DIL 


167.9 




First 1 fb"^ 


Combo 
LJ 
DIL 


171.7 +2:1 

172.2 tli 
166.1 ±5.0 


"•^0 -0.50 
"•>^^ -0.55 




Combo 


171.7 ±2.7 


—0 70 +"'^'' 


Last 0.9 fb"^ 


LJ 
DIL 


170.2 ±11 
175.2 


-0 61 +'''^'' 

U.Di _o.64 



Yan events to the tt mass region, resulting in a systematic 
on ISR-FSR modeling 0. Note that unhke in Ref. [3], 
we coherently shift parameters affecting both ISR and 
FSR, as the uncertainties on the two effects should be 
correlated. We measure the uncertainty due to generator 
choice by comparing MC simulated experiments gener- 
ated with HERWIG and pythia. A systematic on different 
parton distribution functions is obtained by varying the 
independent eigenvectors of the CTEQ6m set [5l|, com- 
paring parton distribution functions with different values 
of Aqcd, and comparing CTEQSl with mrst72 [ssj . 
The gluon fusion fraction uncertainty is calculated by 
reweighting the MC samples to increase the fraction of tt 
events initiated by gluons instead of quarks from the 6% 
in the leading order MC samples used for the measure- 
ment to 20%, which is given as the la upper bound on 
the gluon fusion fraction in (27| . Systematic uncertainties 
due to lepton energy scale are estimated by propagating 
shifts on electron and muon energies within their uncer- 
tainties. Background shape systematic uncertainties are 
obtained by varying the fraction of the different types 
of backgrounds in MC simulated experiments. For lep- 
ton±jets backgrounds, we generate further changes in the 
shapes by varying the used in the calculation of hard 
scattering and shower evolution in the range 
4M^. We also substitute the anti-electron QCD model 
for the non-isolated lepton model in MC simulated ex- 
periments. For dilepton backgrounds changing the shape 
of the Drell-Yan sample according to the difference in the 
missing energy distribution observed in data and simu- 
lation gives one systematic effect. We also shift the fake 
model in ways expected to maximally correlate with the 
reconstructed mass. The systematic uncertainty due to 
limited signal MC statistics is taken as the uncertainty 
on the fit to a constant of the residuals obtained in MC 
experiments (Fig. [6]). We study the effects of limited 
background MC statistics using the bootstrap technique, 
where multiple background MC data ensembles are gen- 
erated. The "pileup" systematic is induced by the possi- 
ble mismodelling of the Minimum Bias events in the MC 
samples. These events are found to deposit more energy 
in jets in simulation than expected from data. It has 
been suggested that color effects may cause a systematic 
bias of order 0.5 GeV/c^ which is not accounted in our 
studies js^j. 

The systematic uncertainties are summarized in Ta- 
ble |Vll The total systematic uncertainty is 1.1 GeV/c^ 
for both the combined and the lepton+jets measurement, 
and 3.8 GeV/c" for the dilepton-only measurement. 



VIII. CONCLUSIONS 

We present the first measurement of the top quark 
mass across multiple decay topologies using a joint like- 
lihood fit. Our procedure includes a full treatment of 
correlations between systematics, and does not assume 
Gaussian likelihoods or symmetric errors in the channels 
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FIG. 11: One-dimensional lepton + jets data with density estimates overlaid using Mtop = 172.0 GeV/c^, Ajes ~ 0.0, and a 
full background model. The expected numbers of events are set to the values from the constrained fit. Shown are the 1-tag 
mj'^'^" (a) and rrijj (b) distributions, and the 2-tag mf^''" (c) and rrijj (d) distributions. 



being combined. In 1.9 fb of data, we measure: 

Mtop = 171.9 ± 1.7 (stat. + JES) 

± 1.1 (other syst.) GeV/c^ 
= 171.9 ±2.0 GeV/c^ 

with cross-checks using events from the Icpton+jets and 
dilepton channels separately: 

Mtop = 171.8 ± 1.9 (stat. + JES) 

± 1.1 (other syst.) GeV/c^ 
= 171.8 ±2.2 GeV/c^ 

(lepton+jets only); 

Mtop = 171.2 tli (stat.) 

±3.8 (syst.) GeV/c^ 
= 171.2 tli GeV/c^ 

(dilepton only). 



This measurement increases our understanding of 
physics in the top quark sector, and contributes to tests 
of the mechanism of electroweak symmetry breaking. In 
addition, the analysis methods and tools described in 
this article will be applicable to other measurements at 
the Tevatron experiments, and soon at CERN's Large 
Hadron Collider. 

The precision of top quark mass measurements at the 
Tevatron is approaching limits set by our understand- 
ing of non-perturbative QCD phenomena. Nevertheless, 
some further improvements are expected as CDF accu- 
mulates a factor of 3-4 times more data during run II 
and as studies of important systematic effects provide 
additional constraints on those uncertainties. 
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